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ABSTRACT 

We have carried out millimeter interferometric observations of the Orion 
Molecular Cloud-2 (OMC-2) FIR 3/4 region at an angular resolution of ~ 3" 
- 7" with the Nobeyama Millimeter Array (NMA) in the Hi3C0+ ( J=l-0), ^^CO 
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(J=l-0), SiO (1^=0, J=2-l), and CS (J=2-l) lines and in the 3.3 mm contin- 
uum emission. Submillimeter single-dish observations of the same region have 
also been performed with Atacama Submillimeter Telescope Experiment (ASTE) 
in the ^^CO (J=3-2) and CH3OH {Jk=Tk-Qk) hues. Our NMA observations 
in the H^^CO^ emission have revealed 0.07 pc-scale dense gas associated with 
FIR 4. The ^^CO (J=3-2,l-0) emission shows high-velocity blue and red shifted 
components at the both north-east and south-west of FIR 3, suggesting a molec- 
ular outflow nearly along the plane of the sky driven by FIR 3. The SiO (t'=0, 
J=2-l) and the submillimeter CH3OII {Jk=7k^Qk) emission, known as shock 
tracers, are detected around the interface between the outflow and the dense gas. 
Furthermore, the ^^CO (J=l-0) emission shows an L-shaped structure in the 
P-V diagram. These results imply presence of the shock due to the interaction 
between the molecular outflow driven by FIR 3 and the dense gas associated 
with FIR 4. Moreover, our high angular-resolution (~ 3") observations of FIR 
4 in the 3.3 mm continuum emission with the NMA have first found that FIR 4 
consists of eleven dusty cores with a size of ~ 1500 - 4000 AU and a mass of ~ 
0.2 - 1.4 Mq. The separation among these cores (~ 5 xlO'^ AU) is on the same 
order of the Jeans length (~ 13 xlO^ AU), suggesting that the fragmentation 
into these cores has been caused by the gravitational instability. The time scale 
of the fragmentation (~ 3.8 xlO^ yr), estimated from the separation divided by 
the sound speed, is similar to the time scale of the interaction between the molec- 
ular outflow and the dense gas (~ 1.4 xlO'' yr). We suggest that the interaction 
between the molecular outflow from FIR 3 and the dense gas at FIR 4 triggered 
the fragmentation into these dusty cores, and hence the next generation of the 
cluster formation in FIR 4. 



Subject headings: ISM: clouds — ISM: OMC-2 / FIR 3, 4, 5 
— ISM: outflows — ISM: molecules — radio lines: ISM 



-stars: formation 



Introduction 



The process of so-called "isolated" star formation, where a sii igle star or binary stars 



form in one molecular cloud core, has been studied intensiyelv (jSaito et al.l Il999l . 12001 



Takakuwa et al.ll2004l : |j0rgensen et al.l 120071 : iTakakuwa et al.ll2007a] ). In this mode, spon- 
taneous gravitational collapse of cloud cores initiates the formation of t he protostar at the 



center, and there is observational evidence of such gravitational collapse (JMvers et al. 



Saito et al. 1996: Mardones et al. 1997: Ohashi et al. 



1997: Momose et al. 1998: Takakuwa et al 



1995 
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2007bl ). Another mode of star formation is cluster formation, where stellar clusters form si- 
multaneously or sequentially in one cl oud core. A vast rn ajority (> 90 %) of stars form 
as members of clusters in our galaxy (JLada fc Ladall2003l ). and the understanding of the 
cluster formation is crucial for more comprehensive understanding of star formation. Previ- 
ous observational studies of cluster fo rmation have proposed that external effects on c loud 
cores, su ch as stellar wind, HII regions (|Williams et al.lll994j . ll995l : lHester fc peschll2005r). su- 
pernova (Kobavashi et al.ll2005l : lYasui et al.ll2006l ). and molecular outflows (ISandell fc Knee 
2OOII : iNakano et al.ll2003l ). are often required to trigger cluster formation. Theoretical stud- 
i es also suggest that such external e ffects are of t en re quired to trigger cluster formation 
( Norman fc Silklll980l ). For example, lYasui et al.l ( 20061 ) have found 52 cluster members in 
"Cloud 2" , which locates in the vicinity of the SNR H I shell, GSH 138-01-94. The age of the 
cluster in Cloud 2 (~ 1.0 x 10*^ yr) is younger than that of GSH 138-01-94 (~ 4.3 x 10^ yr), 
and they suggest that the star formation in Cloud 2 was triggered by the interaction with 
the SNR H I shell. In the 850 /xm dust continuum map of the NGC 1333 region, there are 



several c avity and /or she 



outflows (ISandell fc Knee 



stru ctures observed, which are most likely created by molecular 



20011). One dusty c ore, SK-1, is in one of these shells and is grav- 



itationally bound, and ISandell fc Kned (I2OOII ) suggest that this core is a specific example of 
the outflow-triggered star formation. However, these studies are presumptions mainly based 
on the overall morphology, and detailed studies of cluster formation processes, including 
spatially-resolved observations of both triggers and cloud cores and the investigation of the 
physical processes, have been limited. This is mainly due to the relatively far distance to 
cluster-forming regions and lack of interferometric observations of cluster-forming regions in 
millimeter and submillimeter molecular lines as well as continuum emissions. 

In order to investigate the detailed mechanism of cluster formation, we have initi- 
ated millimeter and submillimeter survey observations of the OMC-2/3 region with the 
Nobeyama Millimeter Array (NMA) and Atacama Submillimeter Telescope Experiment 
(ASTE). The OMC-2/3 region is one of the nearest c l uster-forming regions at a distance 



of 45 pc (IGenzel fc Stutzkil Il989l : Ijohnson et al.lll990l : iLada fc Ladall2003l: iNielbock et al. 



2003). and locates at the north of the Orion A Giant Molecular Cloud (JGatley et al. 



1974 



Sakamoto et al.lll994j : IChini et al.lll997l : iLis et al.lll998l ). Mapping observations of the OMC- 



2/3 region in the 1.3 mm continuum emission have identified 26 dusty c ores (IChini et al. 



19971 : INielbock et al.ll2003l ). and 33 dusty cores in the 850 /xm emission (JLis et al.l Il998l 
^^CO (J=l-0) survey observations of molec ular outflows in OMC-2/3 have found nine out- 
flows at an average separation of ~ 0.2 pc (JAso et al.ll2000l : IWilliams et al.ll2003l ). and our 



recent submillir neter ^^CO (J=3-2) sur vey observations of OMC-2/3 have revealed presence 
of 15 outflows ( Takahashi et al. 20081). Extensive mapping observations of the OMC-2/3 



region in the CS (J=l-0) (iTatematsu et al.l Il993f ) and Hi3co+ (J=l-0) lines (lAso et al. 
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2000l : llkeda et al.ll2007l ) have unveiled physical properties of cloud cores in the cluster-forming 
region, which are quantitatively different from those in isolated star-forming regions. Inter- 
ferometric H^^CO"'"(J=1-0) observations of MMS 7 in the OMC-3 region, one of the dusty 
cores, have revealed the disk- like envelope around M MS 7, which is being dispersed by the 
associated molecular outflow (JTakahashi et al.ll2006l ). 



In this paper, we focus on observations of the OMC-2 FIR 3/4 region (including FIR 
4, an d 5: IChini et al.l 119971) F I R 4 is the strong est 1.3 mm dust-continuum source in 



OMC-2 (jjohnstone fc Ballyl Il999l : IChini et al.l 119971 ). In this region, there are three 3.6 



cm free-free eraissioii sources of VLA 11 (co nsistent with FIR 3) , VLA 12 and VLA 13 
( iReipurth et al.lll999l ). and nine MIR sources (JNielbock et al.ll2003l ). suggesting that FIR 4 
is an active cluster-forming re gion. A molecular outflow driven by FIR 3 h as been found in 
the ^^CO (J=l-0, 3-2) lines Jwilhams et al.ll2003l : Eakahashi et al. Iboosh . which is one of 
the candidate triggers of the cluster formation in FIR 4. Our millimeter interferometric and 
submillimeter single-dish observations of the FIR 4 region should shed light on the detailed 
mechanism of the cluster formation, which will be shown in the present paper. 



2. Observations and Data Reduction 



2.1. NMA observations 



We carried out millimeter interferometric observations of the FIR 4 region in the H^^CO"'" 
( J=l-0; 86.754 GHz), ^^CO (J=l-0; 115.271 GHz), SiO (f=0, J=2-l; 86.847 GHz), and CS 
(J=2-l; 97.980 GHz) lines with the Nobeyama Millimeter Array (NMA), which consists of 
six 10 m antennas, during a period from 2005 May to 2007 January. The data in the H^^CO"'" 
( J=l-0), ^^CO (J=l-0), and SiO {v=0, J=2-l) lines were obtained with the FX correlator, 
which was configured with 1024 channels per baseline and a bandwidth of 32 MHz. For the 
^^CO (1-0) and SiO (2-1) data we made 5-channel binning to increase the signal-to-noise 
ratio (= S/N) of the high-velocity line-wing emission. Thus, the velocity resolution in the 
JJ13Q0+, 12C0, and SiO observations is 0.108, 0.406, and 0.539 km s-\ respectively. The CS 
( J=2-l) d ata were obtained with the digital spectral correlator. Ultra Wide Band Correlator 
(UWBC) dOkumuraet al.lboooh . which has 128 channels and a 1024 MHz bandwidth per 
baseline (1 channel = 24.47 km s~^). The hanning window function was applied to reduce 
side lobes in the sampling. Thus, the effective velocity resolution was widened to 48.94 km 
s~^ (corresponding to 2 channels). The CS (J=2-l) emission is detected over the three 
UWBC channels with the hanning window function, and we integrated over these three 
channels (corresponding to 97.88 km s~^) to make the total integrated intensity map. We 
also obtained continuum data at both the lower (87.090 ± 0.512 GHz) and upper (98.418 
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± 0.512 GHz) sidebands with UWBC, simultaneously with the H^3C0+ (J=l-0) data. To 
obtain a higher S/N in the continuum map, the data of both sidebands were co-added 
(effective observing frequency = 92 GHz = 3.3 mm). Using the AIPS package developed at 
NRAO, we adopted both the uniform and natural UV weighting for the continuum imaging. 
For the imaging of the molecular emissions we adopted the natural weighting. Table [1] and 
[2] summarize the parameters for the line and continuum observations, respectively. Since 
the minimum projected baseline length of the H^^CO"*" (J=l-0), ^^CO (J=l-0), SiO (f=0, 
J=2-l), and CS (J=2-l) observations was 2.9, 3.0, 2.9, and 3.2 kA, our observations were 
estimated to be insensitive to structures more extended than b7" (0.13 pc), 5^ (0.1 2 pc), 
57" (0.13 pc), and 51" (0.11 pc) at the 10 % level, respectively Jwilner fc Welchlll994h . The 
overall uncertainty in the flux calibration was estimated to be ~ 15 %. After the calibrations, 
only the data taken under good weather conditions were adopted in the imaging. 



2.2. ASTE observations 



The submillimeter CH3OH {Jk=7k-Qk\ K=-1, 0, 2; 338.408681, 338.344629, 338.72165 
GHz) data in the FIR 4 region have been taken w i th the ASTE 10 m t e lescop e (Atacama 
Submillimeter Telescope Experiment; lEzawa et al.l (120041 ): iKohno et al.l (120041 )) located at 
the Pampa la Bola (altitude = 4800 m), Chile. The observed frequency range covers the 
CH3OH K-ladder (7/^-6/^, K=-2, -1, 0, 2, and 6) from 338.32 GHz to 338.74 GHz, where three 
of those (7_i-6_i, 7o-6o, 72-62) are detected above 3 a level (= 0.51 K). Observations were 
remotely made from the ASTE operation room at Mitaka, Tokyo on 2006 O ctober 7, 10, and 



11, u sing the network observation system N-C0SM0S3 developed by NAOJ (JKamazaki et al. 



2OO5I ). The half-power beam width of the ASTE telescope is 22" at the observing frequency. 
The typical system noise temperature with the 345 GHz SIS heterodyne receiver in the DSB 
mode was 270 - 440 K, and the typical atmospheric opacity at 220 GHz was ~ 0.05 during our 
observations. The temperature scale was determined by the chopper-wheel method, which 
provides us with the antenna temperature corrected for the atmosphere attenuation. As a 
backend, we used four sets of a 1024-channel auto-correlator, which provide us a frequency 
resolution of 31.12 kHz, corresponding to 0.1 km s~^ at the CH3OH {Jk=7k-6k) frequency. 
The On-The-Fly (OTF) mapping technique was employed to map the 5' x 5' region centered 
on FIR 4, which covers the entire outflow driven by FIR 3 seen in the ^^CO (3-2) emission 
( JTakahashi et al. 1120081 ). The telescope pointing was checked every two hours by 5-points 
scans of the point-like ^^CO ( J=3-2) emission from 0-Cet (aj2ooo= 02^* 19*" 20';80 5j2ooo= 
-02° 08' 40'.'7). The pointing errors were measured to be from 0'.'2 to 1"2 during the observing 
run. We also observed NGC 1333 IRAS 2A and compared the observed CH3OH spectra to 
those obtained with the Caltech Submillimeter Observatory (CSO) telescope, which has the 
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same dish size as that of the ASTE telescope (|j0rgensen et al.ll2004l ). and found that the 
main beam efficiency of the ASTE telescope was ~ 0.4. After subtra cting linear baselines , 
the OTF data were convolved by a Gaussian-tapered Bessel function (JMangum et al.l 120001 ) 
with a FWHM of 14" and were resampled onto a 7" grid. Since the telescope beam is a 
gaussian with a FWHM of 22", the effective FWHM resolution in the restored images is 26". 
The "scanning effect" was minimized by c ombining scans alon g the R .A. and decl. direction, 
using the PLAIT algorithm developed by lEmerson fc Graevd (119881 ) . The typical rms noise 
level in the final image is 0.17 K in T^ at a velocity resolution of 0.5 km s~^. 



3. RESULTS 



3.1. i^cO (J=3-2, 1-0) emission 

Figure [1] shows the distribution of the high- velocity blueshifted and redshifted ^^CO 
(J=3-2) and ^^CO (J=l-0) emission in the FIR 4 region, observed with ASTE and NMA, 
respectively. The ^^CO ( 7=3-2) data were take n by our group as a part of our survey project 
in the OMC-2/3 region jTakahashi et al. Ibood l The ^^CO (J=3-2) line profile toward FIR 
4 is shown in Figure [21 From the peak temperature of the ^^CO (J=3-2) emission, we 
determine that the gas kinetic temperature in the FIR 4 region is ~ 52.5 K. The detected 
velocity range of the ^^CO ( J=3-2) and ^^CO ( J=l-0) emission is -9.5 ~ 29.4 km s~^ and 
3.7 ~ 18.3 km s~^, respectively, while that o f the H^^CQ"^ ( J= l-0) emission, which traces 
the dense gas in FIR 4, is 10.4 ~ 12.5 km s~^ (jlkeda et al.ll2007l ). There are both blueshifted 
and redshifted high-velocity components at the north-east (NE) and south-west (SW) of 
FIR 3. These results imply that the high-velocity components in the ^^CO (J=3-2, 1-0) 
emission trace a bipolar molecular outfiow nearly along the plane of the sky driven by FIR 
3. Furthermore, the high-velocity compo nents of th e ^^CQ emission at the north-east of FIR 
3 are associated with H2 knots found by IYu et al.l (Il997l ). which also supports our outfiow 
interpretation. In addition, the velocity range of the ^^CO (3-2) emission is larger than that 
of the ^^CO (1-0) emission. One of the plausible interpretations of the different velocity range 
is that the submillimeter line traces warmer gas in the vicinity of the jet, where the higher 
velocity is expec ted, due to the differe nt excitation condition between the submillimeter and 
millimeter lines (JRaga fc Cabritlll993l ). 



The size of the SW blue lobe measured from the ^^CO ( J=3-2) image is ~ 0.09 pc, which 
is twice smaller than that of the NE lobe (~ 0.19 pc). At the south-west of FIR 3 both the 
^^CO (J=3-2) and ^^CO (J=l-0) emission exhibit blueshifted and redshifted peaks near 
FIR 4, and that these ASTE and NMA peaks appear to be consistent with each other. On 
the other hand, at the north-east of FIR 3 the prominent blueshifted and redshifted peaks 
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observed with ASTE are not seen in the NMA maps, suggesting that these components 
are resolved out by the interferometer. These results imply that the SW component of the 
outflow driven by FIR 3 toward FIR 4 is more compact than the NE component. From our 
NMA observations, the peak position of the SW red lobe is measured to be aj2ooo=05'^ 35™- 
2?!2, (5j2ooo=-05° 09' 59'.'5, and at the tip of the red lobe there is 3.6 cm free-free emission 
(VLA 12, open blue circle). There are two possibilities on the origin of this free-free emission. 
One is that there is another protostar located at VLA 12, and the other is that this free- free 
emission is originated from the ionized jet driven by FIR 3. Since there is no 3.3-mm dust 
continuum counterpart associated with VLA 12 (see § 3.4), the origin of the 3.6 cm free-free 
emission is probably the ionized jet driven by FIR 3. 

Figure [3] shows a total integrated intensity map in the CS (J=2-l) emission. The peak 
of this emission is (aj2ooo, 5j2ooo) = (05'* 35™ 27:0, -05° 09' 56'.'0), ~ 4" north-east of FIR 
4. This emission extends from FIR 3 to FIR 4, similarly to the high-velocity ^^CO (J=l-0, 
3-2) emission shown in Figure [H Thus, one of the possible interpretations on the origin of 
the CS emission is moderately dense gas in the molecular outflow. Higher velocity-resolution 
observations in the CS emission are required to verify the origin of the CS emission in the 
FIR4 region. 



3.2. Dense gas associated with FIR 4 

Figure H^ shows a total integrated intensity map of the H^^CO"'" ( J=l-0) emission in 
the FIR 4 region, observed with the NMA. There exists a ~ 0.07 pc-scale (~ 30") dense gas 
structure associated with FIR 4. This dense gas s t ructu re has also been identified as AC17 



with the single-dish 45m telescope by lAso et al.l (J2000l ). Hereafter we call this structure 



"FIR 4 clump". The peak position of FIR 4 clump (aj2ooo=05''35'"27!0, 5j2ooo=-05° 10' 
03'.'5), estimated by the 2-dimensional Gaussian fitting to the image, is consistent with the 
position of AC17, and the systemic velocity (=Vsys~ 11.3 km s"^) and the line width (~ 
1.1 km s~^) estimated by the Gaussia n fitting to the H^^CO"*" spectrum at FIR 4 are also 



consistent with the single-dish results ( JAso et al.ll2000l ). Figure Hb shows the distribution of 



the blue- (10.5 km s^^ - 11.2 km s'^) and redshifted (11.3 km s"^ - 12.4 km s"^) Hi3co+ 
(J=l-0) emission in FIR 4. The blueshifted component is distributed from south-east to 
north-west while the redshifted component is distributed from east to west. 

We derived the mass of FIR 4 clump (= Mgas) on the assumption of the Local Ther- 
modynamic Equilibrium (LT E) condition and the o ptically-thin H^^CO^ ( J=l-0) emission. 



using the following equation (JTakahashi et al.ll2006l ) 



M, 



gas 



Skc^mD^ Te^ //i5(J+l)(J + 2) 



IGhuH^B^iXim^CO^ 






KJ- p' 



KAiy, 



where k, h, c, Te^., D, X[H^'^CO^], /i, m, u, J, B and Fj, are the Boltzmann constant, 
the Planck constant, the speed of hght, the excitation temperature, distance, fractional 
abundance of H^^CO"'", dipole moment, the mean mass of the molecule (2.33 a.m.u), fre- 
quency, rotational quantum number, rotational constant, and the flux density, respectively. 
We assume that the excitation temperature of the H^'^CO"'" ( J=l-0) line is equal to 52.5 K, 
which is derived from the peak brightness temperature of the ^^CO ( J=3-2) emission (se e 
§ 3.1). We adopt /i = 4.07 Debye as the H^3C0+ dipole mome nt Jflaese fc Wocxy ll979r). 
The H13C0+ (J=l-0) abundance of 4.5 x 10"" is adopted from kso et all J200oh . and the 



derived r nass of FIR 4 clurnp is ^ 10.5 Mq (If we adopt the abundance of 1.4 x 10 ^° esti- 



mated by iTakahashi et al.l (120061 ) , the mass of FIR 4 clump is 
the momentum (= Pgn s) and the internal gas energy (= E, 
( JTakakuwa et al.ll2003al ) ; 



3.6 Mq). We also estimated 
gas J of FIR 4 clump as follows 



P< 



gas 



M,„,sC, 



gas^eff, 



(2) 



E, 



gas 



-MgasC^ff, 



(3) 



^2 






81n2 



K\ 



m 



mobs 



(4) 



where AVfwhm, Ceff, Tr, and m.obs is the FWHM velocity width of the H^'^CO^ line, 
the effective sound speed, the gas kinetic temperature, and the mass mean molecular of the 
observed molecule. Table [3] summarizes our estimates of these parameters. Due to the effect 
of the missing flux in the interferometric observations (~ 83 % of the total flux is missed), 
these values should be consid ered as a lower limit, and in fact, the estimated mass from the 
single-dish data is ~ 35 Mq JAso et al.lboool l 



3.3. Shock tracer of CH3OH {3^=7 k-Qr) and SiO {v=Q, J=2-l) 

We have detected the submillimeter CH3OH {Jk=7k^Qk', K=-1, 0, 2) lines toward FIR 
4 clump with ASTE. Figure [5] (a) and (b) show the distribution of the blue- and redshifted 
CH3OII (Jii:=7o-6o) emission, and Figure [6] shows the spectrum of the CH3OH (Ji^=7o-6o) 
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line toward FIR 4. The line width of the submiUimeter hne is up to 10 km s~^ at the zeroth 
level, and there is slight indication of the asymmetric line profile, with a more prominent 
blue wing (5.7 km s~^ - 9.9 km s~^) than a red wing (13.1 km s~^ - 14.2 km s~^). The peak 
position of the blue and red lobe is (aj2ooo, Sj2ooo) = {5^ 35"^ 2G'.6, -5° 09' 55'.'4) and {5^ 35"^ 
27^1, -5° 09' 55'.'4), respectively. The distribution of the blueshifted CH3OH emission extends 
along the east-west direction, while the distribution of the redshifted emission extends toward 
the direction of the outflow driven by FIR 3, i.e. along SW-NE. 

We have also detected the SiO (f=0, J=2-l) emission toward FIR 4 clump with the 
NMA, as shown in Figure [5] (c). The peak position of the SiO emission is (aj2ooo5 '^J2ooo) = 
(05'^ 35"^ 2Qm, -05° 09' 53'.'0). The SiO emission is blueshifted, point-like, and locates at ~7" 
north of FIR 4. The peak position of the SiO emission is approximately consistent with the 
peak position of the blueshifted CH3OH emission. 



3.4. 3.3-mm dust continuum emission 

Figure [7^ and Uh show naturally- and uniformly-weighted 3.3-mm continuum images 
in the FIR 4 region observed with the NMA. In these images there appears two main com- 
ponents, one is associated with FIR 3 and the other with FIR 4. The over all distribu- 



tion is consistent with that of the single-dish 1.3 mm continuum observations (jChini et al. 



19971 ). The overall ridge distribution of the 3.3-mm emission associated with FIR 4 also 
resembles that of the H^^CO^ ( J=l-0) emission, although there appears north-eastern ex- 
tension in the 3.3-mm emission. The north-eastern extension in the 3.3 -mm emission may 



trace the "wall" of the molecular outflow, as seen in other observations ( Gueth et al.l 12003 



Moriarty-Schieven et al.ll2006l ). However, the high- velocity CO lobes are unlikely to be the 
cloud-core components traced by the H^^CO^ emission, as discussed in section 3.1. From 
these considerations, we suggest that most of the 3.3-mm emission traces the dense-gas 
component in FIR 4, with a contamination from the outflow component to a certain extent. 

FIR 3 was detected by our 3.3 mm dust-continuum observatio ns as a single source , 



although two MIR s ources MIR 21 and 22, are present in FIR 3 flNielbock et al.l 120031 1. 



Nielbock et al.l (120031 ) have suggested that MIR 21 and 22 are binary stars at a separation 
of ~ 3" (= 1350 AU). However, we consider that MIR 21 and MIR 22 are reflection nebulae 
evacuated by the outflow from FIR 3, because the 3.3 mm peak at FIR 3 locates between 
MIR 21 and MIR 22 and the alignment of MIR 21 and 22 is consistent with the direction of 
the molecular outflow. 

The dusty component associated with FIR 4 is most likely the same identity traced by 
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the H^^CO^ (J=l-0) emission, that is, FIR4 clump. Moreover, the dusty component in 
FIR 4 has subpeaks. These peaks are probably substructures (= cores) in FIR 4 clump. We 
identified those cores with the following criteria; (1) the peak intensity should be higher than 
4 a noise level, and (2) the "valley" among diff'erent cores should be deeper than 1 a. We 
identified five cores in the naturally-weighted image, and six more cores in the uniformly- 
weighted image. Then the total number of the identified cores is eleven. Hereafter, we call 
these smaller-scale cores FIR 4a,b,. . . ,k, as labeled in Figure [7)d. As we mentioned in § 3.1., 
there is no 3.3-mm continuum source toward VLA 12. 

The deconvolved size of these dusty cores was estimated by the 2-dimentional Gaussian 
fitting to the image. The mass of the cores (= M^ust) was derived from the total 3.3 mm 
fiux, Fu, on the assumption that all the 3.3-mm continuum emission arises from dust and 
that the emission is optically thin, using the formula, 

Mdust = — , ., (5) 

KyBi,[Td) 

where we have adopted a value of the mass opacity, Ky = 0.1 f ^^"^"^ j cm^ g~^ 



( jHildebrandl 119831 ) and (3=2. IChini et al.l (119971 ) estimated the value of /3 to be 2 by the 



SED fitting toward FIR 1 and FIR 2 in the OMC-2 region. Since FIR 4 locates on the 
same molecular filament of OMC-2 as FIR 1 and FIR 2, we consider that (3 = 2 is probably 
acceptable for FIR 4. For the dust temperature we adopted Trf=52.5 K, which is derived 
from the ^^CO (J=3-2) emission (see § 3.1). The mean gas density (= n) in these dusty 
cores was derived by assuming a spherically-symmetric shape as follows; 

^= , r-— .3 > (6) 

where Dmaj and Dmin are the deconvolved size along the major and minor axis. The 
typical size, mass and the density of these cores are estimated to be ~ 5 " (= 2250 AU), 0.78 
Mq, and 1.6 x 10^ cm~^, respectively. The total mass of FIR 4 clump estimated from the 
3.3-mm dust continuum emission is 9.4 Mq, which is consistent with that estimated from 
the H^^CO"^ emission (10.5 M©). The average density of FIR 4 clump is 9.8 xlO^ cm~^, 
assuming the spherically symmetric core with a radius of 15" (6750 AU). Table H] summarizes 
these physical properties of the identified cores. 

From the naturally- and uniformly- weighted interferometric 3.3 mm fiux in FIR 4, the 
850 /im fiux is estimated to be 1.2 and 0.58 Jy, respectively, on the assumption of /? = 2, while 
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the 850 //m flux measured with JCMT is 7.5 Jy (jjohnstone fc Ballylll999l ). Then, ~ 84 and 
92 % of the total flux are missed in the naturally- and uniformly-weighted interferometric 
3.3 mm continuum images. Then, the estimated mass of these cores (= M^iust) should be 
considered as lower limits. In Table HI we also list the mass of these cores corrected for the 
effect of the missing flux (= Mcor), which can be considered as upper limits of the mass. 



4. Discussion 



4.1. Interaction between the molecular outflow driven by FIR 3 

and FIR 4 clump 

Several pieces of evidence for the interaction between the molecular outflow driven by 
FIR 3 and FIR 4 clump were found in our observations. The extent of the SW outflow 
lobe from FIR 3 observed in the CO (1-0, 3-2) lines (~ 0.09 pc) is shorter than that of the 
NE lobe (~ 0.19 pc). This suggests that the propagation of the SW outflow is hampered 
by the material at the head of the outflow, and in fact FIR 4 clump locates at the tip of 
the south-western outflow. These pieces of morphological evidence suggest that the south- 
western outflow component is dammed by FIR 4 clump. We also detected the SiO and 
submillimeter CH3OII emission at the tip of the south-western outflow. These molecular 
lines are often o bserved toward shocked regions, caused by t he interaction between outflows 
and dense gas (jAvery &: Chiaol Il996l : iBachiller et al.l I2OOII ). Hence, the detection of the 
SiO and submillimeter CH3OH emission in FIR 4 is probably the chemical evidence for the 
interaction between the molecular outflow and the dense gas. In addition, the line width of 
the H^^CO^ emission at FIR 4 {dYpwHAi ~ 1.1 km s~^) is larger than the average value 



in other OMC-2/3 cores {dYpwHAd ~ 0.8 km s ^; llkeda et al.l 120071 ) which may suggest the 



larger turbulence due to the interaction with the outflow. 

Figure [Hti'-d show Position- Velocity (P-V) diagrams in the ^^CO, H^^CO^, SiO, and the 
CH3OH emission. The P-V diagrams in the ^^CO, 11^^00+, SiO, and the CH3OH emission 
show successive change of the line width along the direction from FIR 3 to FIR 4. First, 
the ^^CO (1-0) emission shows a distinct L-shaped structure in the P-V diagram, and the 
location of the line broadening is at the most upstream. Second, the CH3OH emission also 
shows a similar L-shape with a less broadening, which locates slightly at the downstream 
from the CO L-shape. The SiO emission shows the most extreme blueshifted components, 
which locates at the peak of the blueshifted CH3OH component. Finally, the dense gas 
component traced by the H^^CO^ emission locates at the most downstream of the successive 
emission distribution and is pinched by the L-shaped CO emission distribution in the P-V 
diagram. These successive emission distributions and the velocity structures are likely to 
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trace the detail of the interaction. 

These observational results suggest an interconnection among the molecular outflow, 
dense gas, and the shock traced by the SiO and CH3OH emission. We consider that these 
pieces of evidence support the presence of the interaction between the molecular outflow and 
dense gas. 



4.2. Outflow-triggered star formation 

In the last section, we suggest that there is an interaction between the outflow driven 
by FIR 3 and dense gas of FIR 4 clump. On the other hand, in § 3.4 we have demonstrated 
that FIR 4 clump consists of eleven dusty cores. From these results, we presume that the 
interaction causes the fragmentation of FIR 4 clump into these cores and eventually the 
formation of the next generation of the cluster members. Hereafter we will examine this 
possibility from consideration of the Jeans instability, time scale, and the virial analyses. 



4-2.1. Jeans instability 

First, we will examine whether the fragmentation into the dusty cores can be caused by 
the Jeans instability, by comparing the average separation of the cores to the Jeans length. 
The average 3-din iensional separat ion (= Al) among these cores can be estimated by the 



following formula (jPeng et al.lll998l ): 



Al = 2(— - l)r, (7) 

''•core '^ core /o\ 

7 =-{7 , (8) 

>^ clump 

where 7 , r, ricore, Vciump, and Vcore are the volume filhng factor, the average radius 
of the cores, the number of cores, and the average volume of FIR4 clump and a single 
core, respectively. Vcore was estimated from the average of the projected size of the dusty 
cores (Table H]), and Vciump was estimated from the radius of 15", on the assumption of 
the 3-dimensional spherical symmetry (see eq. 6). The estimated Vcore, Vdump, and r are 
3.1 X 10^^ cm^, 6.5 x 10'^'^ cm^, and 1.1 x 10'^ AU, respectively. Then, the average separation 
Al is estimated to be ~ 5 x 10^ AU. 
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In order to estimate the Jeans length (= Aj), we use the following equation (JNakamura fc Li 



2003); 



2 _ ^Clff 

where Ce//, G, and po are the effective sound speed, gravitational constant, and the 
average density of FIR 4 clump, respectively. Ce// is the effective sound speed in FIR 4 
clump (~ 0.62 km s~^) derived from equation (4), and po is estimated to be 9.8 x 10^ cm~'^ 
in § 3.4. With these values, the Jeans length \j is estimated to be ~ 13 x 10^ AU, which 
is on the same order of the average separation of the dusty cores. Hence it is possible that 
the fragmentation of FIR 4 clump into cores was caused by the Jeans instability. The Jeans 
instability in FIR 4 clump may be triggered by the interaction with the outflow from FIR 
3, and in the next section we will discuss whether the fragmentation was caused after the 
interaction between the outflow driven by FIR 3 and FIR 4 clump. 



4-2.2. Time scale of the fragmentation of FIR 4 clump into cores 

We can estimate the time scale of the fragmentation on the assumption that the frag- 
mentation time scale (= t fragmentation) is the souud crossing time; 

- ^^ 

'^fragmentation ~p^ • l-'-^j 



c, 



eff 



Then r fragmentation IS estimated to be ~ 3.8 x 10^ yr. Next, we estimate the time scale 
of the interaction between the outflow and dense gas, on the assumption that the interaction 
time scale Tinteraction is simi lar to the dynamical tirn e r^ of the north-eastern outflow driven by 



FIR 3, i.e. Tinteraction ~ ^d- iTakahashi et al. I (120081 ) estimated r^ to be 1.4 x 10"^ yr. Then the 
interaction time scale Tinteraction is estimated to be 1.4 x 10^ yr. Therefore, the fragmentation 
time scale (3.8 x 10^ yr) is similar to the interaction time scale (1.4 x 10^ yr), and we 
suggest that the fragmentation was triggered by the interaction with the outflow from FIR 
3. Moreover, the free-fall time in this region (~ 1.3 xlO^ yr), which is estimated from the 
average density of each core (~ 1.6 x 10^ cm~^), is also similar to the interaction time-scale, 
and hence it is possible that Cl ass O/I protostars are already formed in the course of the 
interaction iNielbock et al.l (120031 ). In fact, there exists an MIR source (MIR 24) identified as 



a Class protostar, and it is possible that MIR 24 is formed in the course of the interaction. 
We suggest that in this particular case of FIR 4 the interaction between the outflow 
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and the dense gas most likely triggers the fragmentation into cores, while we cannot exclude 
other possibil ities such as the indu ced fragmentation due to the interaction with the HII 
region (M43') (IHester fc Deschll2005[) as well as spontaneous fragme ntation due to the gravity 
( jlnutsuka &: Miyamal 119971 ) and turbulence (JLi fc Nakamura]l2004l ). 



4-2.3. Outflow-triggered star formation 

Finally, in order to examine whether these dusty cores have a potential to form protostars 
inside, we compare the H2 mass of the cores estimated from the 3.3 mm dust- continuum flux 
and the virial mass. The virial mass (= M^riai) is estimated by the following equation; 



M 



virial 



5rC^ 
2G '■ 



;ii) 



where r, Ce// and G are the radius of the core, the effective sound speed (equation 4) 
and the gravitational constant, respectively. Here, we assume that the H^^CO"'" emission is 
associated with all the dusty cores. Then we adopted the H^'^CO^ line width [^VpwHM] see 
Table H]) at each peak of the dusty core as the line width in the dusty cores. The estimated 
virial masses are within the range of the expected H2 mass of these cores, that is, from M-dust 
to Mcor (Table HI). These results suggest that the observed dusty cores have a potential to 
form stars. 

From the above considerations, we suggest that the interaction between FIR 4 clump 
and the outflow driven by FIR 3 caused the fragmentation of FIR 4 clump into cores and 
triggered the next generation of cluster formation. Our scenario of the star formation in the 
FIR 4 region is summarized in Figure as; 

[STEP 1]: FIR 3 was born and drove the outflow; 

[STEP 2]: The outflow driven by FIR 3 started interacting with FIR 4 clump; 

[STEP 3]: The interaction caused the fragmentation of FIR 4 clump into cores; 

[STEP 4]: These cores will form stars. 

The present stage of the star formation in the FIR 4 region may be between STEP 3 
and STEP 4. We speculate that dusty cores will form protostars inside, and eventually, next 
generation of the stellar cluster in the FIR 4 region. 

P ossible outflow-triggere d star formati on has also been reporte d observationally in NGC 
1333 JSandell fc Kneel boOlh . L1551 NE JVokogawa et all l2003l ). and NGC 2264 IRS 1 
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flNakanoet all 120031 ). In NGC 1333, ISandell fc Kned fl200l[ ) found that one of the cores 
in the 850 /im emission (SK-1) locates at the tip of the shell evacuated by the outflow, and 
they suggest that SK-1 is a possib le example of the outflow-triggered star formation. In 
L 1551 NE, lYokogawa et al.l (120031 ) have found an arc-shape structure of dense gas open 
toward South- West, or the direction of the molecular outflow driven by LI 551 IRS 5, in the 
CS (3-2) emission. They interpret that the dense gas in L1551 NE is affected by the outflow 
from L1551 IRS 5, creating an arc-sh ape, and the protostellar formation of L1551 NE. In 
NGC 2264 IRS l. lNakano et aP J2003h have found that a dense-gas shell seen in the Hi3C0+ 
emission locates around the outflow driven by IRS 1, and at the inner edge of the dense shell 
there are three compact cores. They suggest that these cores are entrained or compressed 
material formed by the outflow and that these cores may represent future sites of the forma- 
tion of a group of the low-mass stars. Our observational studies in the OMC-2 FIR 3/4 region 
have suggested that the interaction between the outflow and dense g as can cause the forma- 
tion of stellar clusters. Recent 3-dimensional MHD simulations by iNakamura fc Lil (120071 ) 
suggest that protostellar outflows can trigger cluster formation through shock compression, 
providing a theoretical support to our scenario. 



5. SUMMARY 

We have carried out high angular-resolution (~ 3" - 7") millimeter interferometric ob- 
servations of OMC-2 FIR 3/4 with the NMA in the Hi3C0+ (J=l-0), ^^CO (J=l-0), 
CS (J=2-l), and SiO {v=0, J=2-l) lines as well as in the 3.3 mm continuum emission, 
and submillimeter single-dish observations with ASTE in the ^^CO (J=3-2) and CH3OH 
{Jk=^k^^k) lines. The main results of our new millimeter and submillimeter observations 
are summarized as follows: 



1. We suggest that the outflow driven by FIR 3 interacts with the 0.07 pc-scale dense gas 
associated with FIR 4 (= FIR 4 clump), from morphological, kinematical, chemical 
and physical evidence. In the morphological evidence, we found that the length of 
the SW outflow lobe from FIR 3 observed in the CO (1-0, 3-2) lines (~ 0.09 pc) 
is shorter than that of the NE lobe (~ 0.19 pc), and FIR 4 clump is located at the 
tip of the SW molecular outflow. These results suggest that the SW lobe has been 
dammed by the interaction with FIR 4 clump. In the chemical and physical evidence, 
we detected the SiO (f =0, J=2-l) emission and the submillimeter CH3OII {Jk=7k- 
6k] K=-1,0,2) emission around FIR 4 clump, which are known to trace shocks caused 
by the interaction between outflows and dense cores. In the kinematical evidence, 
the P-V diagrams in the ^^CO, CH3OH, SiO, and H^^CO^ emission show a successive 
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change of the line width along the direction from FIR 3 to FIR 4. The ^^CO (1-0) 
emission show a distinct L-shaped structure in the P-V diagram, and the location of 
the broadening is at the most upstream. The CII3OH emission also shows a similar 
L-shape at the downstream from the CO L-shape, and the SiO emission shows the 
most blueshifted component at the peak of the blueshifted CH3OH component. The 
dense gas component traced by the H^^CO^ emission locates at the most downstream 
of the successive emission distribution and is pinched by the L-shaped CO emission 
distribution in the P-V diagram. These results indicate the presence of the interaction 
between the outflow and dense gas in FIR 4. 

2. In the 3.3-mm continuum emission, we have first resolved FIR 4 clump into eleven dusty 
cores with a size of 1500 - 4000 AU and a mass of 0.2 - 1.4 Mq. The 3-dementinal 
separation among these cores, 5000 AU, is on the same order of the Jeans length, 
13000 AU. Moreover, the estimated time scale of the fragmentation into the cores, 
3.8 xlO^ yr, is also similar to the time scale of the interaction between FIR 4 clump 
and the outflow from FIR 3, that is, 1.4 xlO^ yr. Therefore, it is possible that the 
fragmentation of FIR 4 clump into the cores has been caused by the Jeans instability 
triggered by the interaction with the outflow driven by FIR 3. We suggest that in this 
particular case of FIR 4 the outflow interaction most likely triggers the fragmentation 
into cores. 

3. We speculate that the dusty cores in FIR 4 will form protostellar sources eventually, 
and that the FIR 4 region is in the course of cluster formation. We suggest that 
the outflow driven by FIR 3 triggered the next generation of the cluster formation in 
the FIR 4 region, and that the FIR 4 region is one of the promising samples of the 
outflow-triggered cluster formation. 
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Table 1. Parameters for the NMA observations 



parameter 



Hl3CO+ (J=l-0) 



^CO (J=l-0) 



SiO (v=0, J=2-l) 



CS {J=2-l) 



Configuration" 

BaseUne [kA] 

Phase reference center {J2000) 

Primary beam HPBW [arcsec] 

Synthesized Beam HPBW [arcsec] 

Velocity resolution [km s"-*^ ] 

Gain calibrator'' 

Bandpass calibrator'^ 

System temperature in DSB [K]'' 

Rms noise level [Jy beam ^^] 



D, C and AB 
2.9-115 

77" 

9^54x5'.' 15 

0.108 km s^i 

3C273, 0420-014, 3C84, 0528+134 

100-400 K 

1.2x10^1 Jy beam -^ 



D 
3.0-30.8 

QJ2000 = 05'^35™26!47,<5 
62" 
6y59x 5'.' 79 
0.406 km s^i 



D 

2.9-23.4 

J2000 = -05°10'00"4 

77" 

10'.'22 X 6^37 

0.539 km s^^ 



0528-fl34 
3C273, 3C84 3C273, 3C84 

150-300 K 100-350 K 

5.7x10^1 Jy beam -^ 2.0x10"^ Jy beam 



D, C and AB 
3.2-113.6 

70" 
6^68 X 4'.'00 
48.94 km s-^ 

3C273, 0420-014, 3C84, 0528+134 

100-400 K 

5.1x10^^ Jy beam ^^ 



^D and AB are the most compact and sparse configurations, respectively. 

'^A gain calibrator, 0528+134, was observed every 20 minutes. 

'^The Bandpass calibrations vifcrc archived by 40 minute observations of each calibrator. 

•^The system noise temperature of the SIS receiver was measured toward the zenith. 



to 

o 
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Table 2. Parameters for the continuum observations 



parameter Fig[7K Fig[33 

Baseline [kA] 2.9-115 2.9-115 

Weighting Natural Uniform 

Beamsize (HPBW) [arcsec] 6'.'96x 4'.'03 (P.A.=-34.59°) 6'.'55x3''34 (P.A.=-40.31°) 
P.A. of the beam [°] -34.59 -40.31 

Rms noise level [Jy bcam^^] 1.2xl0~^ 1.4x10^^ 
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Table 3. Parameters of FIR 4 clump derived from the H^'^CO+ ( J=l-0) emission 



parameter 


Mgas 

[Mq] 


Vsys 

[km s-i] 


AYfwhm 
[km s-i] 


Pgas 

[Mo km s-l] 


[Mq km2 s-2] 


value 


10.5 


11.3 


1.1 


10.5 


5.2 



Table 4. Identified dusty cores 



source 


aj2000 


/3 ,72000 


^maj X ^min 


P.A. 


Ma^st" 


Mcor'* 


n 


^i-virial 


AVpwHM 


Ce// 








[" X "] 


["] 


[Mq] 


[Mq] 


[xlO^ cm-3] 


[Mq] 


[km s-i] 


[km s-i] 


FIRS 


05'' 35"" 27!6 


-05° 09' 34'.'0 


11.7x3.6 


140.46 


3.0 


18.7 


3.1 








FIR 4a" 


05 35 27.4 


-05 09 49.0 


5.4x3.5 


34.05 


0.51 


3.2 


1.8 


1.1 


0.53 


0.47 


FIR 4b" 


05 35 26.6 


-05 09 51.0 


6.3x3.9 


150.20 


1.1 


6.9 


2.6 


2.6 


1.28 


0.68 


FIR 4c 


05 35 27.0 


-05 09 54.0 


9.0x4.0 


168.88 


0.45 


5.9 


0.59 


1.6 


0.57 


0.48 


FIR4d 


05 35 26.2 


-05 09 55.9 


4.7x4.4 


139.87 


0.76 


9.6 


2.2 


1.2 


0.55 


0.47 


FIR4e 


05 35 26.9 


-05 09 57.5 


8.4x5.7 


4.302 


1.0 


14.3 


0.93 


3.8 


1.31 


0.69 


FIR4f 


05 35 26.6 


-05 09 58.9 


6.1x3.1 


168.26 
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^Identified in the naturally-weighted map. 

''The FWHM size estimated by the 2-dimensional Gaussian fitting to the image. 

'^Derived from the total fiux density integrated within the 3cr contour levels, with D=450 pc, /3=2, and T^=52.5 K. See texts for details. 

•^112 mass corrected for the effect of the missing flux. See texts for details. 

•^Assuming a spherically-symmetric core with a radius of 15" (=6750 AU). 



to 

00 



-24- 



ASTE 



o 



Q 10' 00 



o 

Q 10' 00 



11' 00 



35ni3is 



12C0 (7=1-0) 
3. 7~6. 1 km/s 




29s 27s 25s 
RA (J2000) 



29s 27s 25s 
RA(J2000) 



Fig. 1. — Distribution of the high- velocity blueshifted (blue contour) and redshifted (red 
contour) ^^CO (J=3-2) (left) and ^^CO (J=l-0) (right) emission in the 0MC-2/FIR 4 
region. Panel (a) and (c) show the distribution of the blue (-9.5 km s~^ - 4.6 km s~^) and 
red (16.4 km s"^ - 29.4 km s^^) lobe in the ^^CO (J=3-2) line observed with ASTE (?), while 
panel (b) and (d) show the distribution of the blue (3.7 km s~^ to 6.1 km s~^) and red (15.1 
km s~^ to 18.3 km s~^) lobe in the ^^CO ( J=l-0) line observed with the NMA, respectively. 
Open black circle s show the field of vie w of the NMA observations, green crosses the position 
of MIR sources (INielbock et al.ll2003l ). blue circles the position of 3.6 cm free-free sources 



S 



( iReipurth et al.l Il999fl. and red squares the peak position of the 1.3 mm dust- continuum 



emission ( Chini et al.l 119971 ). Filled ellipses at the bottom left corner of each panel indicate 
the beam size. Contour levels of these maps start at ± 5 a levels with an interval of 5 a. 
The rms noise levels (1 a) of the panel (a), (b), (c) and (d) are 1.2 K, 0.45 Jy, 1.2 K, and 
0.70 Jy beam~^, respectively. 
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Fig. 2.— ^^CO (J= 3-2: 345.705 GHz) spectrum at FIR 4 taken with ASTE 



(JTakahashi et al. 1 120081 ) The rins noise level (la) is 0.5 K in T^ at a velocity resolution 
of 1.1 km s^^. A dashed line shows the systemic velocity (~ 11.3 km s~^), obtained from 
the Gaussian fit to the spectrum. 
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Fig. 3. — Total integrated intensity map of the CS ( J=2-l) emission in the 0MC-2/FIR 4 
region, obtained with the NMA. Symbols in the Figure are the same as in Figure [H Contour 
levels of this map start at ± 3 o" levels with an interval of 1 a. The rms noise level (1 a) of 
this map is 5.1x10"^ Jy beam~^. 
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Fig. 4. — Total integrated intensity map of the H^'^CO^ (J=l-0) emission (left), and maps of 
the blueshifted (10.5 - 11.2 km s"^) the redshifted (11.3 - 12.4 km s"^) H^^CO+ components 
(right) in the FIR 4 region, obtained with NMA. Symbols in the Figure are the same as in 
Figure [H Contour levels of these maps start at ± 3 a levels with an interval of 1 a. The 
rms noise levels (1 a) in panel (a), the blue and red contour in panel (b) are 0.81, 0.49 and 
0.49 Jy beam~^, respectively. 
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Fig. 5.— ASTE CH3OH ( J^=7o-6o) and NMA SiO (t;=0, J=2-l) maps in the FIR 4 region. 
Panel (a) and (b) show the distribution of the blue (5.7 km s~^ - 9.9 km s~^) and red (13.1 
km s~^ - 14.2 km s~^) lobe in the CH3OH (Jr-=7o-6o) line, respectively. Panel (c) shows 
the distribution of the blue lobe in the SiO (1^=0, J=2-l) line at a velocity range from 3.4 
km s~^ to 6.1 km s~^. Symbols in the Figure are the same as in Figured! Contour levels in 
panel (a) and (b) start at ± 5 a levels with an interval of 5 a. Contour levels in panel (c) 
start at ± 3 cr levels with an interval of 1 a. The rms noise levels (1 cr) in panel (a), (b) and 
(c) are 0.024, 0.015 K, and 0.11 Jy beam~^, respectively. 
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Fig. 6.— CH3OH ( J^=7o-6o; 338.408 GHz) spectrum at FIR 4 taken with ASTE. The rms 
noise level (la) is 0.17 K in T^ at a velocity resolution of 0.5 km s^^. A dashed line shows 
the systemic velocity (~ 11.3 km s^^) in the FIR 4 region estimated from the NMA H^^CO^ 
spectrum. 
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Fig. 7. — 3.3 mm dust-continuum maps in the FIR 4 region observed with the NMA. The 
left panel shows the naturally-weighted map, and the right panel the uniformly-weighted 
map. Symbols in the figure are the same as in Figure [H Contour levels start at ± 3 a noise 
levels with an interval of 1 a. The rms noise level (1 cr) in panel (a) and (b) are 1.2x10"^ 
and 1.4x10"^ Jy beam~^, respectively. 
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Fig. 8.— Position-Velocity (P-V) diagrams of the ^^CO (J=l-0) (upper left), Hi3C0+ 
( J=l-0) (upper right), SiO {v=0, J=2-l) (lower left) and CH3OH (J/^=7o-6o) (lower right) 
emission in the FIR 4 region. Two horizontal lines in each panel show the position of FIR 
3 and FIR 4, while a vertical line the systemic velocity of 11.3 km s~^ in FIR 4 estimated 
from the H^^CO"^ data. The cut line of panel (a), (b), and (d) is adopted to be along P. A. = 
30°, passing through the position of FIR 3 and FIR 4. The cut line of panel (c) is taken to 
be along P. A. = 35°, passing through the compact SiO emission and FIR 3. Contour levels 
in panel (a) and (d) start at ± 5 a levels with an interval of 5 a. Contour levels in panel 
(b) start at ± 3 cr levels with an interval of 3 a. Contour levels in panel (c) start at ± 2 cr 
levels with an interval of 1 a. The rms noise level (1 a) in panel (a)-(d) are 0.57 Jy beam~^, 
0.12 Jy beam~^, 0.2 Jy beam"^, and 0.17K, respectively. 
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Fig. 9. — Schematic picture of the outflow-triggered star formation in the FIR 4 region. In 
STEP 1, FIR 3 was born and drove the outflow. In STEP 2, the outflow driven by FIR 3 
interacted with FIR 4 clump. In STEP 3, the interaction between the outflow and FIR 4 
clump caused the fragmentation of FIR 4 clump into cores. In the final STEP, these cores 
form stars. The current stage of the FIR 4 region may be between STEP 3 and STEP 4. In 
these figures, horizontal lines show the plane of the sky and the observers are located at the 
bottom. 



